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We show that QCD4 with transverse confinement can be approximately compactified into QCD2 
with a transverse quark mass that is obtained by solving a set of coupled transverse eigen- 
value equations. In the limits of a strong coupling and a large number of flavors, QCD2 further 
admits Schwinger QED2-type bosonized solutions. We therefore examine phenomenologically the 
space-time dynamics of produced particles in string fragmentation by studying the Wigner function 
of produced bosons in Schwinger QED2, which mimics many features of string fragmentation in 
quantum chromodynamics. We find that particles with momenta in different regions of the rapidity 
plateau are produced at the initial moment of string fragmentation as a quark pulls away from an 
antiquark at high energies, in contrast to classical pictures of string fragmentation with longitudinal 
^ , space-momentum-time ordering. 
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I ^' I. INTRODUCTION 

Soft particle production by string fragmentation is an important process in high energy collisions. In the strongly- 
O • coupled regime in nucleus-nucleus collisions at RHIC energies, the dynamics of string fragmentation falls within the 
1^ ] realm of non-perturbative quantum chromodynamics (QCD), and the detailed mechanism of the production process 
. cannot yet be described from the first principles of QCD. Phenomenological particle production models based on 
preconfinement [ll, parton-hadron duality cluster fragmentation string-fragmentation Q , dual-partons , 
, the Venus model W, the RQMD model [7|, multiple collision model [1], parton cascade model color-glass 
J> ' condensate model . the AMPT model 1^, the Lexus model [l^, and many other models have been put forth 
Q^ I and are successful in describing some aspects of the data. These models incorporate some features of QCD, but their 
fundamental foundation in non-perturbative QCD is still lacking. It is of interest to study the physics of particle 
2^ _ production in string fragmentation with QCD-inspired. low-dimensional non-perturbative models. 

. . Two-dimensional quantum electrodynamics [ij, [l^ [H, [13, [H, [H, [13, [U, 113, HI, HI] in one space and one time 
' coordinates (QED2) furnishes an interesting arena for such a study. This is a quantum mechanical system in which 
a fermion interacts with an antifermion through a linear gauge field potential in the Coulo mb gaug e. It is a system 
where the Higgs phenomenon occurs and the global chiral symmetry is spontaneously broken [14| Il5l [l6l. Il7|. When a 
positive fermion and a negative antifermion are separated in such a system, the vacuum is so polarized that the positive 
and the negative charges are completely screened, in a manner similar to the confinement of quarks. The quantum 
field theory of massless fermions interacting with a gauge field in QED2 can be solved exactly Q M, M, [13, H [H ■ 
The interacting field turns out to be equivalent to the quantum field theory of free bosons with a mass. It was 
cd ' demonstrated by Casher, Kogut and Susskind that when a fermion separates from an antifermion in the limit 
of infinite energies in QED2, the rapidity distribution of the produced particles will exhibit the property of boost 
invariance. For a finite energy system, the boost-invariant solution turns naturally into a rapidity plateau, whose 
width increases with energy as In(Y^) [23l. l2^. 

Rapidity plateau distributions of produced particles have indeed been observed in high-energy e'^-e~ annihilation 
which produces initially a separating quark and an antiquark [25, 26, 27, 28, 29, 30]. Rapidity plateau distributions of 
produced particles have also been observed in pp collisions at ^/s ~ 200 GeV by the BRAHMS Collaboration [IH . The 
close similarities of the features of empirical fragmentation data and theoretical QED2 fragmentation results motivate 
us to examine in this paper the circumstances in which the longitudinal dynamics of a transversely-confined QCD4 
system can be approximately compactified into QCD2. We also need to examine further how multi-flavor QCD2 can 
be approximated by QED2. 

Our renewed interest in the non-perturbative production mechanism arises from the puzzling observations of the 
momentum distribution of particles associated with a near-side jet. In high-energy heavy-ion collisions, a near-side jet 
is characterized by the presence of associated particles within a narrow cone along the trigger particle direction. Its 
characteristics resemble those of a jet in pp and peripheral heavy-ion collisions. In addition to this 'jet component', 
there is an additional 'ridge component' of associated particles at A(j)^0 with a ridge structure in A77, where Acj) and 
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A77 are the azimuthal angle and pseudorapidity differences measured relative to the trigger particle [H, [H, [13, [H, 

M, M, m, M M M, m, SI, M m M m M M MMM- 

In the phenomenon of the ridge associated with the near-side jet, it is observed that (i) the ridge particle yield 
increases with the number of participants, (ii) the ridge yield appears to be nearly independent of the trigger jet 
properties, (iii) the baryon to meson ratios of the ridge particles are more similar to those of the bulk matter than 
those of the jet, and (iv) the slope parameter of the transverse distribution of ridge particles is intermediate between 
those of the jet and the bulk matter [33|-[13- These features suggest that the ridge particles are medium partons, at 
an earlier stage of the medium evolution during the passage of the jet. The azimuthal correlation of the ridge particle 
with the jet and the presence of strong screening suggest that the associated ridge particle and the trigger jet are 
related by a collision. A momentum kick model was put forth to explain the ridge phenomenon [s^ . |54, 55 , 5a, [57| . 
The model assumes that a near-side jet occurs near the surface, collides with medium partons, loses energy along its 
way, and fragments into the trigger and its associated fragmentation products (the "jet component"). Those medium 
partons that collide with the jet acquire a longitudinal momentum kick along the jet direction. They subsequently 
materialize by parton-hadron duality as ridge particles in the "ridge component" . They carry direct information on 
the momentum distribution of the medium partons at the moment of jet-(medium parton) collisions. Based on such 
a description, the experimental ridge data of the STAR, PHOBOS, and PHENIX Collaborations [H, [H, H S, SI 
reveal an early parton momentum distribution that has a thermal-like transverse momentum distribution and a 
rapidity plateau extending to rapidities as large as |?;| ^4. 

The early stage of a nucleus- nucleus collision consists of simultaneous particle production processes involving quarks 
pulling away from antiquarks (or qq diquarks) at high energies. As experimental data of elementary production 
processes have a rapidity plateau structure [1^, [H, [23, [H, [H, [sO, [3l| , partons at the early stage of the nucleus- 
nucleus collision can possess the characteristics of a rapidity plateau. A rapidity plateau would also be expected 
from the theoretical model of string fragmentation in QED2 Jisl. [Tol . [23l . [23 |. the dual-parton model [^, the Bjorken 
hydrodynamics [l9l |. the Lund model 01, the Venus model [a|^color-glass condensate model the Lexus model 
[13|, and many other models. Thus, the possibility of a rapidity plateau for the ridge particles should not come 
as a surprise. However, the fact that a parton of large absolute rapidity can occur together with a jet of central 
rapidities, as revealed by the PHOBOS data and the momentum kick model analysis, poses an interesting conceptual 
question. For those medium partons with large magnitude of the longitudinal momentum to collide with the jet so as 
to become an associated ridge particle in the PHOBOS experiment, the partons must be present in the longitudinal 
neighborhood of the jet at the moments of the jet-medium collision, at an early stage of the nucleus- nucleus collision. 

In many classical descriptions of particle production processes, a particle with a large absolute rapidity is associated 
with a large separation from the longitudinal origin. For example, in the Lund model of string fragmentation 
there is a momentum-space-time ordering of produced particles in the fragmentation of a string. Consider the case in 
the center-of-mass system. Those particles with the smallest absolute rapidity will be produced nearest to the origin 
at the earliest times, while those with the largest absolute rapidity will be produced farthermost from the origin at 
the latest times [3, [^ . 

In Bjorken hydrodynamics [l^, the spatial configurations are invariant with respect to a boost in rapidity. There is 
a longitudinal momentum-space ordering of the produced particles. Particles at the local longitudinal location z ~ 
have rapidity zero and it is necessary to go to a relatively large absolute longitudinal coordinate |z| to find another 
particle with a large absolute rapidity \y\. 

In another classical model of Landau hydrodynamics [H, [s^, [6^, one starts with a fluid initially at rest. As the 
fluid evolves, there is also an ordering of the rapidity with the longitudinal coordinate and time. It is necessary to go 
to larger absolute longitudinal coordinates at later times to find particles with larger absolute rapidities. 

In these classical pictures of produced particles, there is a space-momentum-time ordering. Those partons with 
large values of absolute rapidities are not produced at the early stage of the j et p roduction and therefore cannot 
be associated with a jet by a collision. On the other hand, the PHOBOS data [43 and the momentum kick model 
analysis [ssl. [s^. [ssl. [Sa. [STj indicate that a central rapidity jet and another parton at a large absolute rapidity can be 
associated by a collision. How does one resolve the difference of these two seemingly contradicting results? 

It should be pointed out that the momentum-space-time ordering of produced particles in the above classical 
descriptions is obtained in flnal-state-interaction models. There is however another class of well-justified initial-state- 
interaction models, in which the constituent particles or p roduced particles are present at or before the collision process 
at t = 0. Notable examples are the parton model [6ll | and the Drell-Yan process [G^I, which assume the presence 
of constituent partons inside the hadron or in the quark- antiquark sea before collisions. The multi-peripheral model 
[63t assumes the exchange of a tower of rcggeons and the plateau of particles are present before the collision at t=0. 
The dipole approach in photon-induced reactions [s^l assumes the fiuctuation of the incident photon into a dipole 
quark-antiquark pair before its collision with the hadron. 

With regard specifically to the dynamics of produced particles in the fragmentation of a string, it is important to 
point out that there are important quantum effects [l^ . [20l . [65j arising from the vacuum structure in strongly-coupled 
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non-pcrturbativc quantum field theories that are beyond the realm of classical considerations. The proper space- 
time dynamics should be based on a quantum field theory where the quantum effects of the vacuum structure play 
appropriate roles. As is well known, in theories where the fermions are weakly interacting and can be isolated, the bare 
vacuum suffices for the description of the dynamics in a weak-coupling perturbation theory. In contrast, however, in 
our systems with strong coupling between fermions, such as in massless QED2 and similar field theories, the vacuum 
of the filled Dirac sea is the appropriate vacuum which participates in the dynamics. The self-consistent dynamics of 
the filled Dirac sea in the non-perturbative response to the presence of gauge field and fermion disturbances lead to 
peculiar quantum effects of the axial anomaly [l^, [13, HO, [IHl , the confinement of the fermion , and the appearance 
of massive bosons [l^ . The particle production process in string fragmentation should therefore be treated in a 
non-perturbative quantum field theory framework. 

While a complete non-perturbative theory based on full quantum chromodynamics is desirable, it is however not 
available. Simplifying approximations are needed. We shall show that a system of transversely-confined QCD4 can 
be approximately compactified to a system of QCD2 with a transverse quark mass that is obtained from a set 
of coupled transverse eigenvalue equations. We shall show that in the limit of strong coupling and large number 
of flavors, a multi-flavor massless QCD2 admits Schwinger QED2-type solutions with an effective coupling constant 
that depends on the number of flavors. In the absence of rigorous non-perturbative QCD4 solutions, we are therefore 
justified to examine phenomenologically the dynamics of the system in the solvable quantum field theory of QED2, 
on accounts of its close theoretical connection with the transversely-confined QCD4 and its desirable properties of the 
proper high-energy rapidity plateau behavior, confinement, charge screening, and chiral symmetry breaking. 

The space-time dynamics of produced particles in string fragmentation can be obtained by evaluating the Wigner 
function of the produced particles when a strongly-coupled fermion separates from an antifermion at high energies 
in QED2. In contrast to the classical description of particle production, we shall find that produced particles with 
different momenta in different regions of the rapidity plateau are present at the moment when the overlapping fermion- 
antifermion pair begin to separate. 

The space-time dynamics of partons in the central rapidity region has been obtained in the color-glass condensate 
model In contrast to the strongly-coupled non-perturbative field theory discussed here, the color-glass condensate 
model considers the region of a; ^ 1 and pr ^ Aqcd in the weak-coupling limit of isolated quarks and gluons. The 
space-time dynamics has been evaluated in the boost-invariant limit of infinite energies with the valence quark in 
the no-recoil approximation. The description is appropriate for the production of particles in the central rapidity 
region with large transverse momenta. In the present case we wish to examine the space-time dynamics of the 
production of particles with large (pseudo)rapidities (|Ar/| ^ 4) and small transverse momenta {pt > 0.035 GeV), 
such as those ridge particles detected by the PHOBOS Collaboration [4^, in a system of finite energy. These particle 
are in the region of large x and small pt ^ Aqcd, and the dynamics is within the realm of strong coupling. It is 
therefore appropriate to investigate here the space-time dynamics of produced particles in string fragmentation in a 
strongly-coupled non-perturbative quantum field theory. 

This paper is organized as follows. In Section II, we spell out explicitly the circumstances in which a transversely- 
confined QCD4 can be approximately compactified to a massive QCD2 of fermions with a transverse mass tot which 
is the eigenvalue of transverse confinement. In Section III, we study how the multi-fiavor massless QCD2 and QED2 
are related in the limit of a strong coupling and a large number of flavors Nf. The dynamics can then be considered 
as those of QED2 with an effective coupling constant that is proportional to ^JNf. In Section IV, we review the 
particle production processes in massless QED2 and show how the Wigner function can be obtained from the initial 
fermion current. In Section V, we give an explicit example in the fragmentation of a string with a finite energy and 
compare the experimental rapidity distribution data in high energy e^-e~ annihilation with the theoretical rapidity 
distribution. In Section VI, The Wigner function of the produced particles is evaluated as a function of time to study 
the space-time dynamics of produced particles in the fragmentation of a string in high energy e"*"-e~ annihilation. In 
Section VII, we generalize the results for a single string to the case of many identical strings at different locations and 
show how the Wigner function can be evaluated. In Section VIII, we give numerical examples of the Wigner function 
in the fragmentation of identical strings. The Wigner function exhibits the effects of interference between different 
strings. In Section IX, we discuss schematically the space-time dynamics in the case of the fragmentation from many 
independent string in the context of high-energy heavy-ion collisions. We present our discussions and conclusions in 
Section X. 



II. APPROXIMATE COMPACTIFICATION OF TRANSVERSELY-CONFINED QCD4 TO QCD2 

Under certain appropriate conditions, quantum chromodynamics in 4 dimensions (QCD4) can be approximately 
compactified to quantum chromodynamics in 2 dimensions (QCD2). It is useful to describe the circumstances and 
assumptions leading to this simplification. Previously, t'Hooft showed that in the limit of large Nc with fixed g^Nc 
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in single-flavor QCD4. planar diagrams with quarks at the edges dominate, whereas diagrams with the topology of a 
fermion loop or a wormhole are associated with suppressing factors of \/Nc and respectively [66j . In this case 

a simple-minded perturbation expansion with respect to the coupling constant g cannot describe the spectrum, while 
the l/iVc expansion may be a reasonable expansion, in spite of the fact that Nc is equal to 3 and is not very big. The 
dominance of the planar diagram allows one to consider QCD in one space and one time dimensions (QCD2) and the 
jhysics resemble those of the dual string or a flux tube, with the physical spectrum of a straight Regge trajectory 



67 1 . Since the pioneering work of t'Hooft, the properties of QCD2 systems have been investigated by many workers 
M HI, m, [ei,!?^, [nl, UA III, Izi, lli, Izi, III, H Izi, iO, im. The flux tube picture of the longitudinal dynamics is 



further phenomenologically supported by hadron spectroscopy |82| and the limiting average transverse momentum of 
produced particles in high-energy hadron collisions (ssj . 

The description of the QCD4 dynamics as a flux tube and the dominance of the longitudinal motion over the 
the transverse motion in string fragmentation provide the motivation to compactify, at the least approximately, the 
transversely-confined QCD4 to QCD2 and QED2. Such a link was presented briefly in [sHi- We would like to reiterate 
the derivation further to show how the coupling constant of a transversely-confined QCD4 can be related to the 
coupling constant in QCD2 in such a compactification. 

We start with QCD4 in which the dynamical variables are the quarks fields ■0°'' and the gauge fields where 
a and b are the color and flavor indices respectively, and v = 0,1,2,3. The color and flavor indices will not be 
displayed explicitly except when they are needed. The gauge fields are coupled to the quark fields which are in turn 
coupled to the gauge fields, resulting in a coupled system of great complexities. The transverse confinement of the 
flux tube can be represented by quarks moving in a nonperturbative scalar field m(r), part of which may contain the 
flavor-dependent quark rest masses. The equation of motion of a quark field ip is 

{ilp-m{r)}i' (1) 

where 

i^=7''n^, (2) 

n^^PM + s^A*' (3) 

A, = r'^Al, (4) 

and r° are the generators of the color SU(iVc) group. The equation of motion for the gauge field is 

D^F^^^ ^g^l^Yri^. A^, - 0, 1, 2, 3, (5) 

where 

Ft,y = df,A^ - dyAf, - ig[Af„ A^], (6) 

D^F^'' ^d^F'^'' -ig[A^,F^''], (7) 

F,,=t'^F^,. (8) 

The fermion current that is the source of the gauge field is 

f = giJl^Tip. (9) 

We shall consider the problem of string fragmentation in which the momentum scale for the longitudinal dynamical 
motion of the source current is much greater than the momentum scale for the transverse motion, \v'^\ 3> \v^\, l^^l, 
where v is a typical source velocity. In the Lorentz gauge, the associated gauge field A^^ is proportional to (l,v). 
Under the dominance of the longitudinal motion over the transverse motion in string fragmentation, A^ and A^ can be 
neglected in comparison with the magnitudes of A^ or A'^ . For the flux tube configuration, it is further reasonable to 
assume that the gauge fields A^ and A'^ inside the tube depend only weakly on the transverse coordinates r = (x^, a;^). 
It is meaningful to investigate these fields inside the tube by averaging them over the transverse profile of the flux 
tube. After such an averaging, A^ and A^ can be considered as a function of (xP^x'^) only. As a consequence, the 
equation of motion for the quarks become 

{7^00 + 7^1 + 1-P2 + 7'n3 - m(r) } ^ = 0, (10) 
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where Ho and have been approximated to be independent of the transverse coordinates. To separate out the 
longitudinal and transverse degrees of freedom, we expand the quark field in terms of spinors /z^ [23j 



iPix) = (GlMl - G2M2)/+ + iGl^l3 + G2M4)/- 



(11) 



where Gi_2 are functions of the transverse coordinate r, f± are functions of (x", x^, and the spinors have been chosen 
to be eigenstates of the gamma matrix a with = Mi, 2 and ayU3,4 = — [23[, 



Ml 



M2 




M3 




and 



M4 



Working out the Dirac matrices, we obtain the set of coupled equations of motion for the quarks 

- m(r)Gi/+ + (-pi + ip^)G2f+ + (H^ + n^)Gi/_ 
-(pi + *p2)Gi/+ + m(r)G2/+ + (n" + n-^)G2/- 
(n" - n3)Gi/+ - TO(r)Gi/_ - (pi - ip2)G2/- 
-(n" - n3)G2/+ + (pi + *p')Gi/_ - m(r)G2/- 

Applying 11° — on (I13a[) and using (|13cp and (jl3dp , we get 

[(n" - n3)(n" + n^) - m^r) - p|]Gi/_ = [(-pi + zp2)m(r)]G2/- 



0, 
0, 
0, 
0. 



(12) 



(13a) 
(13b) 
(13c) 
(13d) 



(14) 



By the method of the separation of variables, we can introduce the eigenvalue such that the above equation can be 
separated into an equation for the transverse coordinates and an equation for the longitudinal and time coordinates. 



[pI + m2(r) - m|]Gi(r) = [(pi - ip2)m(r)]G2(r), 



(15) 



[(n" - n3)(n" + if) - n4]f^{x",x^) = o. 



(16) 



By applying 11° — if on (jl3bp and using (|13c[) and (|13dp . we can separate out the same longitudinal equation but 
the following transverse equation. 



[p| + m2(r) - TO|]G2(r) = -[(pi + ip^)m{r)]Gi{r). 



(17) 



Similarly, by applying II" + 11^ on (|13cp and using (|13ap and (jl3b[) . we can separate again the same transverse equation 
(|15p for Gi , and the following longitudinal equation 



[(n" + n-^)(n" - n^) - to|]/+(.t",x3) = o. 



(18) 



The set of four coupled equations of (|13ap - (jl3d|) from Eq. ([T]) (or PTI)) ) for the quark therefore becomes the set of 
four equations of (fT5l) - p8)) . Eqs. (fT5|) and pT)) provide the equations for the solution of the eigenvalue under 
the boundary condition that the functions Gi and G2 are confined, with vanishing probabilities at |r| — > 00. Note 
that Eqs. psp and (jl7p is applicable not only to an azimuthally symmetric transverse potential but also to an 
azimuthally-asymmetric potential m{r). For an azimuthally symmetric transverse potential m{r), we can further 
write the transverse wave functions as 



Gi(r,0) = 9i{r)e^-t 
G2{rA) = 52(r)e'^(''+i)^. 
Eqs. and (fTTl) then become a coupled set of equations for gi{r) and (72 (j') [s^ - 

(JTTL { V] 

[V>l{y)+nf{r)-ml]gi{r) = -^^^52(0, 
[p2(i/) +m2(r) -m^] £,2(0 = i^$^3i(0> 



(9r 



where 



P^(-) 



1 i9 i9 J/ 

r c)r dr r 



(19a) 
(19b) 



(20a) 
(20b) 

(21) 
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In solving the transverse eigenvalue equations (|15p and (fTT)) . (or (|20ap and ([20b|) with the boundary condition 
of transverse confinement, an eigenstate is characterized by transverse quantum numbers and the corresponding 
eigenvalue tot depends on these quantum numbers. Thus, the transverse quark mass will depend on the transverse 
quantum state of the quark in the flux tube. Excitation of the quark in the flux tube will lead to a greater transverse 
mass. There is thus an additional degree of freedom for the transverse quark mass. The excitation of the transverse 
state will lead to a broader width of the transverse momentum distribution for the quark. Composite objects formed 
by pairing a quark with an antiquark will likewise acquire a broader transverse momentum distribution [83| . 

Eqs. (|16|) and are the longitudinal equations for a quark in the two-dimensional space-time a;'^). They can 
be rewritten as a Dirac equation in the two-dimensional space-time as follows. We introduce a two-dimensional quark 
spinor 

V'2D = UA (22) 



and use two-dimensional gamma matrices [7J, |81| 

72D = h ) ' = h ) ' 72D72D = 72D = 1 g -1 ) ' ^ ' 

Then Eqs. and can be written as 

{{Uy (tf)^ +7|D72"D[^^^"] - to|} = 0. (24) 

This same equation ([M)) can also be obtained from 

{7°Dno + T^oHg - tot} V'2D = {72d(?'o + gAa) + 7|d(P3 + .9^3) - wt} V'2D = 0, (25) 

which is the equation of motion for a quark in two-dimensional gauge fields of Aq and A3. 

Note that the above equation depends on .94Dj4^(4d) where we have added the subscript '4D' to g and to indicate 
that the coupling constant is the coupling constant in QCD4, and the gauge fields are from the Maxwell equation 
([5]) in the four-dimensional space-time. In the two-dimensional space-time of QCD2, the corresponding equations of 
motion for the quark and the gauge fields are 

{72d(P0 + g2DAo(2D)) + 72d(P3 + ^20^3(20)) " ™t} '020 = 0, (26) 

and 

-C'/x^ITd = 320-020720 'rV'2D = ^20^ = 0, 3, (27) 

where are given by Eqs. ([6l)-(l8|) with A,^(4d) replaced by Aj^(20) and by 520- 
To cast Eq. ()25|) in the form of (|26p we can identify 

540^^(40) = 320 Ap(20)- (28) 

To cast Eq. ([5]) in the form of Eq. (P7)) . we need a relationship between and jjj^^. With the quark field as given 
by ini) and ([22]), we have 

fiD = .940 ■04O 740 '^■040 

= .g4o(|Gi(r)|2 + |G2(r)nv;2o72D^02O, 

= (.94o/32o)(|Gi(r)p + |G2(r)n.j-j,. (29) 

This indicates that the j'^^i current depends also on the transverse spatial density distribution. As our focus is in 
the longitudinal dynamics for high-energy string fragmentation, it suffices to average the quark transverse density 
|Gi(r)p -I- |G2(r)p over the flux tube transverse profile to relate j40 approximately with j2o as 

(j4o). = (54o/.92o)((|Gi(r)|2 + \G2{v)\y ^f^,,, (30) 

where the transverse averaging of O is defined as 

(O), = / dvO{\G^{v)\^ + \G2{v)\^). (31) 
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The source term of the MaxweU equation involves a coupUng constant and the gauge field quantity in the Dirac 
equation also involves the coupling constant. The equations of motion of the quarks and gluons (|26|) and ^ can be 
cast in the forms of ([26]) and ((27l) by transversely averaging Eq. ([5]) over the profile of the flux tube and by relating 
the coupling constants by the following renormalization 

52D-54D((|Gi(r)|2 + |G2(r)|2))^. (32) 

We can get an estimate of the relation between 520 and by considering the case of a uniform transverse flux tube 
profile with a transverse radius Rt, 

(|Gi(r)p + |G2(r)p) ^ e(i?T - \r\)/nR^T- (33) 

We have then 

32D = 91d/{^Rt) = ^c^s/Rl- (34) 

If we consider the case of e+-e^ annihilation at ^/s = 29 GeV, the average p\ of produced pions is 0.255 GeV^ [l^ 
corresponding to a flux tube radius of order 

Rt ^ - 0.4 fm. (35) 

A typical case of a = 0.4 and a flux tube radius of Rt = 0.4 fm then gives 

glj,^OAOGcY\ (36) 
On the other hand, our earlier comparison relates with the string tension coefficient b in the linear potential is 

52D - 26, (37) 

which also give ~ 0.4 GeV^, for a string tension of 6 = 0.2 GeV^ = 1 GeV/fm. The estimate of the QCD2 
coupling constant given above in Eq. p4|) is therefore consistent with the string tension coefficient. This relationship 
also gives a relation between the string tension b and other physical quantities, 

b = 2ajRl. (38) 

Thus, with the many approximations outlined above, a transversely-confined QGD4 can be compactified as massive 
QCD2 with a transverse quark mass tot obtained by solving the transverse eigenvalue equations (|15p and (jl7p . 

Our explicit formulation to compactify a four-dimensional space-time in a flux tube to a two-dimensional space- 
time serves to provide a great simplifications which allows us to examine the dynamics of systems dominated by the 
dynamics along the longitudinal direction, as in high energy string fragmentation. The transverse degrees of freedom is 
subsumed by the use of the transverse quark mass rriT , and the longitudinal and the transverse degrees of freedom are 
decoupled. If needed subsequently, the transverse properties may be approximately restored by using the transverse 
eigenfunctions and their corresponding transverse momentum distributions of the quarks and antiquarks. For example, 
while we get the rapidity distribution of the produced bosons, its transverse distribution may be considered by treating 
the boson as a composite object made up of a quark and an antiquark, for each of which the transverse momentum 
distribution can be obtained from the transverse eigenfunctions Gi(r) and G2(r), as discussed previously in [sst . 



III. RELATION BETWEEN QCD2 AND QED2 



In the last section we have shown that the transversely-confined QCD4 can be approximately compactified to QCD2 
with a transverse quark mass determined by the transverse confinement of the quarks within the flux tube. We can 
examine the flavor degrees of freedom and write down the flavor index b in the QCD2 equations of motion (|26p and 
(j27p . In the general case when the transverse mass depends on the flavor (as the current quark masses may be 
flavor-dependent), the QCD2 equation of motion for the quark is 

{7°(Po + 52D^o)+7'(P3 + 320^3) -m^}^A' = 0, 6=1, ...,7V/, (39) 

and the Maxwell equation for the gauge field is 



6=1 



(40) 
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where the subscripts '2D' of 7, ^, Ay, and F^^'^ have been henceforth omitted, for brevity of notation. 

We ideahze to the case of a system with flavor symmetry so that the transverse masses and the quark wave functions 
of the underlying system are independent of the flavor. Then the gauge field equation (j40p becomes 

j:,^pt^i^ = Nfg2Ti^^YT^^, (no sum in h). (41) 

As the strength of the source of the gauge field increases with Nfg2i), the case of a large flavor number TV/ and a 
large coupling constant is expected to be in the realm of strong coupling. In this case, the multi-flavor QCD2 can be 
best investigated by bosonization as it display a theory of a scalar particle with weak self-interactions. The bosonized 
massive QCD2 gives rise to a sine-Gordon equation with periodic symmetry that does not yield simple solutions. In 
this case of strong coupling and in high energy processes for which the energy of the system is much greater than the 
transverse masses, it is reasonable to consider massive QCD2 in the massless limit and treat the mass as a perturbation 
in the mass-perturbation theory [13, HI] . 

Accordingly, we can consider the presence of a perturbation of the gauge field A'^"' {v = 0, 3) in massless QCD2 
and can choose the perturbation to be along the direction of a color component a. The corresponding field strength 
tensor F^"^" will also be along the a direction and commutes with A'^"-. The equation of motion for the gauge field 
(|4T|) becomes the Abelian Maxwell equation 



df, [d'^A'"' - d^'A'"'] = Nff^"''. (42) 

where 

j"'"'' g2DV'S''-r''V'^ (no sum in b). (43) 

Upon taking into account gauge invariance and the singularity of the Greens function, the induced current j'^''^'' for 
this Abelian case of massless quarks is related to the gauge field perturbation by [3, 0| 



■i^.ab _ 9 2D 



J ' = 

The above Abelian Maxwell equation is satisfied if 



A''- -d-'^d^A^^^ 



(44) 



- DA""^ - ^ll^A""- = 0. (45) 

TT 

The equation of motion for the gauge field perturbation then becomes the Klein-Gordon equation for a boson with a 
mass square given by 

= Nfglj,/^. (46) 

as obtained previously b y F rischmann, Sonnenschcin, Trittmann and their collaborators, and confirmed by numerical 
calculations [tH. FtgI. ItsI. l85l| . These authors argue further that considerations such as given above for the classical 
equations of motion in massless QCD2 is applicable especially in the large Nf limit, for which the partition function 
is dominated by the classical configuration in the functional integral. Investigation in this limit captures the quantum 
behavior of the theory. 

For each color index a of there is a massive Klein-Gordon equation (|45|) . The massive bosons are therefore 
[N'l — l)-fold degenerate, which is the degeneracy of the gluon. The bosons are the quanta of the excitation arising 
from the perturbation the color gauge fields A^"", they have the color octet property of the underlying gauge field 
perturbations. Furthermore, the coupling constant (^20 is proportional to as given by Eq. (j34p . The boson mass 

M as given by Eq. is proportional to 1/ . A boson mass of this magnitude is also expected for a gluon in the 
four-dimensional space-time confined to a flux tube of radius i?^ . Based on the above, the massive boson states are 
effectively gluons because (i) it has the proper degeneracy of the gluon, (ii) arises from the excitation of the color 
gauge fields A^°- , and (iii) it is of the right order of magnitude of the gluon mass as in QCD4 space-time in a confined 
flux tube of radius . 

One can keep track of the degrees of freedom from a different perspective as one compactifies QCD4 into QCD2 
and QED2. In the asymptotic freedom regime of QCD4, light quarks and gluons are free and essentially massless at 
short distances. However, for large distance behavior in the transversely-confined QCD4 in a flux tube, gluons and 
quarks are confined in the tube and they acquire a mass of order H/Rt corresponding to the zero-point energy in the 
confining tube. When transversely-confined QCD4 is compactified into QCD2, the degree of freedom are the massive 
quarks and the gauge fields A° and A"^, where one of the gauge fields can be eliminated by fixing the gauge and the 
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other can be written in terms of the fermion currents. It may appear at first that the gluon degrees of freedom have 
disappeared. However, as we note above, the excitation of the gauge fields A'^"' leads to massive color-octet boson 
states with the properties of confined and massive gluons. Thus, the gluons in QCD4 reappear as color bosons in the 
excitation of the gauge field. These massive bosons can also be considered as excitations accompanied by the coupling 
of quark-antiquark pairs in the color-octet configuration. The excitation of the A'^"' gauge fields producing one of 
these bosons leaves the remaining fermions in the corresponding color-octet configuration such that the combined 
system of the gluons and quarks is color-neutral. 

What is the fate of these color bosons (or gluons) in the picture of string fragmentation and jet-parton interaction 
as examined in [H, [13, HI, [H, 113'^ One can envisage the picture that as they emerge and decay into the detected 
colorless particles, the color of these color scalar bosons can be evaporated with the emission of soft gluons that carry 
the color away from one boson to another neighboring boson until all bosons and remaining fermion composite objects 
become colorless. Such a color evaporation picture then leads to the model of parton-hadron duality 0] that have 
been extensively used in the discussion of high-energy collision processes as for example in the Glassma model [llj . 

The bosonization of multi-flavor QCD2 in the larg e Nf limit indicates that the QCD2 can be approximated as 
QED2 with three essential modifications [7l|, [tI, [tsTFtqI. Firstly, in the limit of massless QCD2, the boson mass 
square is given by 320^//"" with an additional factor of Nf, and this massive boson state possesses color and has the 
degeneracy of {N^ ~ !)• Secondly, the non-zero transverse quark mass in massive multi- flavor QCD2 leads to a mass 
term which gives rise to interactions between the color bosons. In addition to modifying the particle spectrum [86j . 
the interactions arising from the non-zero transverse quark mass will affect the dynamics of particles after production. 
Thirdly, the magnitude of the fermion current generating the excitation of the gauge field increases linearly with the 
flavor number Nf. As a consequence, the number of produced boson particles should also increase linearly with the 
flavor number Nf. 

While the spectrum is known in the limit of large Nf and Nc for massless QCD2 [tiI, [tI, [t^ [t^ , the spectrum for 
the realistic case of TV = 3 and Nf ~ 3 in massive multi-flavor QCD2 has not been obtained. We shall carry out a 
phenomenological fragmentation study by following the Abelian bosonization of QCD2 to QED2 and shall equate the 
mass m of the boson in QED2 as the experimental transverse mass of detected pions in our analysis. We shall also 
likewise treat the QCD fragmentation problem as a single flavor QED2 problem with the multiplicity of produced 
particles modified by the multiplicative factor of Nf. 

We note in passing that we can carry out an additional coupling constant renormalization for a multi-flavor QCD2 
with flavor symmetry by identifying 

dm = 92D\/n]- (47) 
We introduce rescaled gauge field A'^ related to by 

9mA'^= 92DAf,. (48) 

Then Eqs. and (gl]) becomes 

{7°(Po+.92D^o)+7'(P3+.92D^3)-™t}^' = 0; 6 = l,...,iV/, (49) 

and the Maxwell equation for the gauge field becomes 

D^iF'Y" = g'2iy4'^YTi}^ (no sum in h). (50) 

The above set of equations and ([50]) is the same as the set of single-flavor QCD2 equations with an effective 
coupling constant g'^^,. The dynamics of the flavor-symmetric multi-flavor QCD2 is therefore the same as a single- flavor 
QCD2 with a renormalized coupling constant 17213 given by g2D\/Nf in (j47p . 

IV. WIGNER FUNCTION OF PRODUCED PARTICLES 

After discussing the relationship between QCD4 and QED2, we shall briefly review and summarize previous QED2 
results in order to show how the Wigner function of produced particles in string fragmentation can be evaluated. 
In conformity with the usual notation, we denote the coupling constant in QED2 by e which will subsequently 
take on values related to the QCD2 coupling constant — (72D (with the choice of the negative sign by convention) 
or phenomenologically as times the transverse mass of a detected pion. For convenience, we shall label the 
two-dimensional coordinate and momentum hy x = {x°,x^) = {t,x) and p = {p°,p^), respectively. 

From the work of Schwinger [T^ , it is known that QED2 involving massless fermions with electromagnetic interac- 
tions is equivalent to a free boson field (j) with a mass m = e/y/n, where e is the coupling constant. We can understand 
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this remarkable property of massless QED2 from the foUowing hncs of reasoning. (For a simple pedagogical discussion 
of QED2, see Chapter 6 of [II]-) One starts with charged fermions occupying the negative-energy Dirac sea. If there 
is a current disturbance j^^ in some region of space, it will lead to an electromagnetic gauge field A'^. This electro- 
magnetic field will affect all wave functions of all particles, including those in the Dirac sea. The resultant wave 
functions generate a charge current j^. This charge current in turn generates a gauge field A^ which interacts with 
the fermions to generate the current. The coupling of the gauge field A^^ and the fermion current j'' in a self-consistent 
manner is a problem of great complexity. Remarkably, when the gauge invariance property of the current is properly 
taken into account, the generated current jf^ is related to the the introduced local electromagnetic disturbance. This 
current in turn leads to an electromagnetic field A^, through the Maxwell equation. If this generated electromag- 
netic field is self-consistently the same as the electromagnetic field which was first introduced, then we find that 
OA^ is proportional to A^. The electromagnetic field satisfies the Klein-Gordon equation appropriate for a boson 
with a mass, OA^^ + {e'^ /tt)A^^ = 0. In the Lorentz gauge, we can represent the gauge field A'^ in terms of the boson 
field (j) by A^ — e^'^d^(f)/m. QED2 involving interacting massless fermions is therefore equivalent to a field of free (p 
boson field with quanta of mass m. 



The relation between the fermion and boson field quantities is [14 



r 



(51) 



where j'' is the fermion current which can be taken to be a real quantity, and e''" is the antisymmetric tensor 
e — — £01 — — -L- 

In QED2, the string fragmentation process is described as the time evolution of a fermion separating from an 
antifermion at a center-of-mass energy ^/s, with the production of the quanta of the boson field (or the perturbation 
of the A'-'- field) at subsequent times. It can be formulated as an initial-value problem. If the current j'' arising from 
the fermions is initially known, then the subsequent dynamics of the boson field (f> can be inferred at all times. 

The stipulation that the interacting fermion field is exactly equivalent to a free massive boson field facilitates the 
solution of such an initial-value problem. We can represent the dynamics of the system through the evolution of a 
real massive pseudoscalar field </>, with the most general solution 



(x,t) = 



1 



/27r 



(52) 



where p ■ x ~ p^t — p^x. The coefheient c{p^) and its complex conjugation c*(p^) can be obtained from the initial 
fermion currents j^{x, t), and j^{x, t) BXt = 0. By using Eq. (|5ip . we therefore have 



—^d^(j){x, t)]t=o = 0), 



(53) 



and 



From these initial conditions and Eq. (|52p . we obtain 



dt(l){x,t)]t=o^3\x,0). 



(54) 



P 



(55) 



and 



c*(V) 



^P{p')^j'{p') 
P 



(56) 



where j^^{p^) is the Fourier transform of j''(a;,0), 

j7^(pi) = 



dxe-'P ^j''(x,0). 



With the coefficients c(p^) and c*{p^) thus determined, the energy P° of the system can be determined from 



(57) 



P0 = i 
2 



/ 5(/)(x, t)\ / d4>{x, t) 



\ dt 



dx 



{4>{x,t)y 



dx. 



(58) 
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It is useful to separate 4>{x, t) in terms of positive and negative frequency components, 
where 



and 



The energy of the system is then 



27r2p" 
dpi 



c(p )e 



P° = y da; 
This leads to the energy 



dt(b^+^ [x, t)dt4>^~^ [x, t) + (a;, t)d^ct)(-\x, t) + m20(+) (x, t)</)^-) (x, t) 



pO 



dpic(pi)c*(pi)/2, 



(59) 



(60) 



(61) 



(62) 



(63) 



which is clearly a time-independent quantity. The momentum distribution of the bosons is given by [22 

dN c(pi)c*(pi) 



dpi 



2p" 



2pOe2 



V(p')+J^(p') 



P 



where = a/ (p^)^ + ^ri^. The rapidity distribution of the produced particles is [23 



dA^ _ TT 

"dy ~ 2^ 



(64) 



(65) 



We have thus obtained a simple relation between the rapidity distribution and the Fourier transforms of the initial 
fcrmion charge currents. 

To obtain the Wigncr function of the produced particles, we rewrite Eq. as 



P" = / dxdk 



2tt 



dt^^+^ {x + J)at0(-) {x~^) + a.^w {x + J)a.0(-) {^~^) + m^^w {x + ^)0(-) {x - ^) 



(66) 



The total energy of the system is related to the Wigner function f{x,k) by 

P° ^ j dx dk k°f{x,k), 



where = \/W + rri? , and the total number N is 

N = j dx dkf{x,k). 

By comparing Eqs. (j66p with (|67p . the Wigner function of produced bosons is given by 

f{x,k)-- 



ds e^^'^ 



27rfcO 



Upon substituting the explicit forms of 0'-^-' in terms of c{p^) and c*{p^) from Eqs. (pO]) and (|6ip . we obtain 



/(a;,fc) 



47rfc0 



2pV 



(67) 



(68) 



.(69) 



(70) 
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where p° = ■\/(fc + (7/2)2 _|_ ^2 g^j^(j pO' _ {k — + ni^ . We can alternatively rewrite the above as 



rf9C0s{Q[x--(p-'(g)-p"(g))i]}il^^ ^___L^c(fc + |)c*(fc _ |). (71) 

Therefore, if an initial fermion current is given, its spatial Fourier components j^{p^) can then be evaluated using 
Eq. ([57|) . and the spatial Fourier coefficients c{p^) and c*{p^) of the bosons can be obtained from Eqs. ([SS]) and (j56p . 
These Fourier coefficients can be used to evaluate the Wigner function in Eq. ([7T|) . 

As the longitudinal momentum is often represented by the rapidity variable y, with = mcoshy and |k| = A; = 
TOsinhy, we can represent the Wigner function as a function of the longitudinal coordinate x, and the rapidity y as 



/(x,2/) = fc°/(2:,fc) 



(72) 

k—m sinh y 



V. STRING FRAGMENTATION AT A FINITE ENERGY 

As a simple example, we shall apply the above results to particle production in the high-energy annihilation of e+ 
with e^, which falls within the domain of multi-flavor QCD4. The approximate consideration of the dynamics as 
occurring within a flux tube in the fragmentation of a string, allows one to approximate the transversely-confined 
QCD4 to QCD2. Multi-flavor massless QCD2 in turn becomes a set of — 1 Abelian QED2-type systems in the 
large Nf limit, with the production of the Schwinger QED2-type states. Because the boson mass for the realistic case 
of massive multi-flavor QCD2 with Nf ~ 3 is not known, we shall content ourselves with using phenomenologically 
the experimental observed transverse mass of the produced pions as the mass of the produced bosons. To take into 
account the linear dependence of the multiplicity on the number of flavors, we can consider the fragmentation of a 
single-flavor string, with the flavor number ascribed to be the effective charge v — Nf of the separating charges. 

We therefore start with a positive charge with charge ve separating from a negative charge —veaXa finite energy 
-yi in the center-of-mass system, to simulate approximately the dynamics as a quark separating from an antiquark 
after an electron-positron annihilation. At t = the fermion and the anti-fermion charges superimposed so that the 
total charge density of the system is zero: 

/(a:,0) = 0. (73) 

However, the initial charge current j^(x,0) is non-zero, and it is a symmetric function of x. In the case of an infinite 
energy, the current can be represented as a sum of two delta functions of the two sources moving with the speed of light 
out to -l-a; and —x directions (isj and the corresponding terms in the boson field (f> are step functions. This suggests 
a simple generalization of the initial (5- function current for finite energies, in which one replaces their step-functions 
in 0{x ±t) by a smooth function of diffusivity cr psj . 

, Jt:v f 'X + t ^ X — 1\ 

4){x,t) = ~- tanh tanh . (74) 

2 \ a / 

The diffusivity a is related to the total invariant mass -/s of the system given in Eq. ((8T|) below. From Eq. (|5T|). the 
field (/)(a;, t) leads to the current j''(a;, t) ~ 0: 

/(x,t) = ^( i— I ■ (75) 

^ ' ' 2a \^cosh2(^) cosh^(2^) ) 

and 

'"^'''^ = (cosh^W) ^ cosh^W)) ■ ^''^ 

Strictly speaking, the fermion currents should travel with a speed slightly less than the speed of light. At high 
energies, the difference between their speed and the speed of light is so small that it can be neglected. At t = 0, we 
have j°(x, 0) = and 

ji(x,0) = i^— . (77) 

cr cosh (x/cr) 
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For this initial current distribution j'^(x,0) the Fourier transforms of .7^(2:, 0) are 

fip') = 0, (78) 

and 



VTT ^TT V27rsmh(7rpitT/2) 

From Eqs. these initial charge current Fourier components leads to the Fourier coefficient 

sinh(p^7r(T/2) 

We can use Eq. to calculate the energy P'^ = ^/s which is independent of time, we obtain a relation between a 
and ^/s^. 

o^^- (81) 



The rapidity distribution is therefore [25 



where 



2f2 



sinh^ ^ ' 



(82) 



J^^TT^msinhy 

4 = TT-r ■ v83j 

Thus, the rapidity distribution is boost-invariant, dN /dy = , in the limit of very high energies. At a finite energy ^/s, 
the rapidity distribution therefore turns into a plateau structure, with a half-width yo at half maximum characterized 
by ^0 = 1-491 and 

= sinh-i ( . (84) 



The half-width yo varies with ^/s of the fragmenting string as yo esc sinh ^{y/s) ~ In(y^), as expected. 



VI. EVALUATION OF THE WIGNER FUNCTION FOR A SEPARATING FERMION AND 

ANTIFERMION PAIR 

It is of interest to study a specific numerical example to illustrate the space-time dynamics of string fragmentation. 
We can consider the experimental rapidity distribution of charged pions in e~^-e~ annihilation at ^/s = 29 GeV [25l.l26l|. 
In such an annihilation, heavier particles such as kaons and baryons are also produced and they carry a fraction of the 
initial energy. By integrating out the experimental dN/dyptdpt data in Table 11 of Ref. [11], the particle multiplicities 
of the produced pions, kaons, and baryons are in the ratio of iV^ : Nfc : N]\[ = 80 : 15 : 5 and the total energies of the 
produced particles are 19.11, 7.10, and 2.53 GeV for pions, kaons, and nucleons respectively psj . 

Rather than attempting to model the production of heavier mesons here, we shall instead simply consider pions 
only and study an idealized situation in which the system fragments only into pions. We can consider the pions to 
take up =19.11 GeV of the total energy and fit the rapidity distribution of the pions with Eqs. ((82|) and ([83]) 
which contain only two unknown parameters: the mass m and the effective charge v. The pions have an average 
transverse momentum of 0.48 GeV [25| . which corresponds to a pion transverse mass of 171^^^, = \/ pip -V fn\ — 0.50 
GeV. In approximately compactifying QCD4 to QCD2 and QED2, we have subsumed the details of the transverse 
degree of freedom by using the transverse mass of quarks. With the quark mass modified to be the transverse quark 
mass by this compactification, the corresponding produced boson mass should also be modified to be the transverse 
boson mass of the observed boson. Accordingly, the boson mass m of the produced particle in Eqs. ([83]) should be 
taken to be the pion transverse mass of m = m^^ = 0.50 GeV. The theoretical rapidity distribution of charged pions 
can then be calculated with Eq. ([5^ with the effective charge v as the only unknown parameter. The effective charge 
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FIG. 1: (Color online) The rapidity distribution of pions in hi g h-ener g y e^-e annihilation at ^/s =29 GeV, out of which 
■^5^=19. 11 GeV goes into producing pions. The data is from [251. [26l. W\ and the solid curve is from Eq. (|82|) with m — 0.50 
GeV, y = 1.73. 



value oi V ^ [l^ gives a good description of the experimental rapidity plateau of the produced pions [l^, [1^, [131 , 
as shown in Fig. 1. It is interesting that as Nj is related to v by N f ^ . It appears phenomenologically that the 
effective number of flavors participating in the excitation of the vacuum for e^-e~ annihilation at this energy is 3. 
This effective charge is in line with the discussions in in Section III, where the strengths of the underlying color source 
generating the excitation of the gauge fields is proportional linearly to the number of flavors N f. 

From Eq. (I81|l . the width parameter of the initial charge current a is 0.33/GeV or 0.065 fm, which is a narrow 
current distribution. With the knowledge of a and the Fourier coefficient c(p^) can be evaluated the Wigner 
function from Eq. (|71|) by direct numerical integration. Fig. 2 shows the Wigner function /(x, y) as a function of x, 
for various values of y at i =0 in Fig. 2(a), for positive y at t =0.4 fm/c in Fig. 2(b), and for negative y at t =0.4 




-0.6 -0.4 -0.2 0.2 0.4 0.6 

X (fm) 

FIG. 2: (Color online) Wigner function f{x,y) as a function of x for different values of positive y. Fig. 2(a) gives f{x,y) for 
t — 0, and the dotted curve is the approximate analytical result for f{x,y = 0) at t = given by Eq. (|86p . Fig. 2(b) gives 
f{x,y) for positive y at t — 0.4 fm/c, and Fig. 2(c) for negative y a,t t = 0.4 fm/c. 
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fm/c in Fig. 2(c). The range of rapidities in Fig. (1) span over difTcrent regions of the rapidity plateau as one can see 
in Fig. 1. One notes the following interesting features. 

1. At t=0, the peak of all Wigner functions occur at a; = 0, indicating that all produced bosons with different 
momenta in different regions of the rapidity plateau are present at t = 0. This picture is in contrast to the 
momentum-space-time ordering of the classical picture of particle production where particles of larger momenta 
\y\ are produced at larger |a;| at a later time. 

2. The Wigner function is symmetrical with respect to the changes of both x to —x and y to —y. That is, 
f{x,y;t) = f{-x,-y;t). 

3. The Wigner function for the zero-mode with y = behaves approximately as jy^g-smla-l ^^ggg (jg^ below) 
and is positive definite. The width of the dominant Wigner function peak for |y| ^ decreases as a function of 
increasing |y|, as expected from the uncertainty principle. 

4. Except for Wigner function of the zero-mode with y — which is positive definite, the Wigner function of other 
longitudinal momenta oscillates as a function of space and time and can assume negative values at locations 
away from the dominant peaks. 

5. Except for Wigner function of the zero-mode with y = which does not depend on time, the dominant peaks of 
the Wigner function moves to the positive longitudinal direction for positive y and to the negative longitudinal 
direction for negative y. The speed of the movement of the peak position increases with an increasing magnitude 
of the longitudinal momentum. 

From the above features, one can readily understand the movement of the peak positions of the Wigner functions for 
different y values, as they follow the motion of what is expected from classical physics. 

How do we understand the relatively large width of the spatial distribution of the Wigner function of the zero-mode 
with y = 0? The half width of the distribution is about 0.2 fm, which is large compared to the width parameter a ^ 
0.065 fm for the charge current. What physical quantity determines the scale of this width? It is instructive to check 
the Wigner function for this zero mode case. For y ~ 0, the Wigner function becomes 

f{x,y = Q) = - I dgcos(gx) ,„^^^^ . 2.^ 



2n Jo ' ' (g2/4) to2 sinh^(g7rcr/4) 

2 /-oo 



2 /"Oo 2 
V 1 Tfl 

'^^^^^('^-) (,V4)+..^ Q("-'-^')' ^''^ 
where f2 ^ 2^o/'^'^ = 2 x 1.491/7r(T. In the spatial region where fix ^ 1, we have 



f{x,y = 0) / 



- mz/2e-2"l^l. (86) 

We plot the function mi^^e"^"*'^' as the dotted curve in Fig. 2(a). We find that except in the region ot x < 0.05 fm, 
it agrees with the exact Wigner function f{x, y = 0). It is also independent of i, as — p''' = for /c = in Eq. (|7ip . 
The analytical result shows that the width of Wigner function for the zero mode with y = is governed by the mass 
boson mass as (width)~ hjm. 

How do we understand the origin of the large range of momentum y values that are present in the Wigner function 
at t = 0? From our derivation of the Wigner function, we note that the presence of an initial charge current with 
a very narrow width in j^^{x,t) leads to a wave packet of the boson field (f> in Eq. ([52|) . as a coherent sum of a 
large number of waves of different momentum components p^. The greater the energy of the fragmenting string, the 
narrower is the initial spatial current j^{x, t), and the greater will be the spread of the momentum components of the 
boson wave packet (f>{x,t), as indicated by the large width of |c(p^)|. The coherent sum of the boson field gives rise 
to the corresponding Wigner function that depends on the correlation between the Fourier components at different 
momenta. The Wigner function has a spatial extension much greater than that of the initial charge current. 

In physical terms, there are two peculiar quantum effects that make QED2 different from those from classical 
considerations of string fragmentation. Firstly, the produced bosons occur initially in a spatial region more extended 
than the width of the initial fermion current. They are governed more by the produced boson mass and the uncertainty 
principle than by the spatial width of the initial current. Secondly, produced particles of different momenta emerge 
at the initial moment of string fragmentation. These peculiarities does not violate causality because the dynamics 
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of string fragmentation in QED2 is not just a two-body problem involving only the fragmenting valence fcrmion and 
antifermion. It involves a many-body problem of particles, antiparticles, and gauge fields interacting self-consistcntly 
with particles in the Dirac sea over an extended region of space and time, giving rise to the excitations of the vacuum 
which manifests themselves as quanta of the massive boson field [3, [13, Ej] ■ 



VII. INTERFERENCE IN THE FRAGMENTATION OF MANY IDENTICAL STRINGS 

The process of string fragmentation is a quantum phenomenon. We therefore expects effects of quantum interference 
when many identical strings fragment. As a well-defined interference phenomenon, it possesses its own intrinsic 
theoretical interest. Furthermore, the presence or absence of the characteristics of interference may be used to provide 
useful information on the fragmenting strings concerning their properties of being identical or non-identical, when 
they occur in close vicinity of each others. 

We consider the occurrence of n identical strings at longitudinal locations Xi centered around the longitudinal origin, 
in the center-of-mass frame. For simplicity, we shall take the fragmentation to be occurring at the same time and 
that the strings are identical in their characteristics. In this case, the fragmentation can be idealized as n interacting 
fermion-antifermion pairs. In each pair, the fermion and the antifermion pulls apart in opposite directions. The 
time-like fermion current j^{x,0) = is initially zero at all spatial points and the initial space-like current j^(x,0) 
pulls each fermion-antifermion pair away from each other. 

In representing the sum of the fermion currents in a many-string system, we would like to introduce the concept of 
the sense of a fragmenting string. Consider the fragmentation of a string starting with a fermion and an antifermion 
locating at a spatial point. The sense of a string specifies the direction of motion of the fermion, with the antifermion 
traveling in the opposite direction. We choose the convention that the sense C is equal to 1 when the fermion travels 
to the positive longitudinal direction and C is -1 when the fermion travels to the negative longitudinal direction. 
The sense of motion clearly is immaterial for a single string. However, they arc important in the fragmentation of a 
many-string system, as we shall see. 

The total current j^^^(x,t) from the n fermion-antifermion pairs can be described as 



1 



1 



2cr ^ \ cosh' 

2 — 1 \ 



2 (x-Xi+t^ 



cosh' 



2 fx--2a^^ 



(87) 



and 



1 , , ev v-^ 



V-^hM^) cosh^(^: 



(88) 



where the width parameter a will need to be self-consistently determined in terms of the total energy -y/s of the 
n-string system because the strings will interfere with each other in generating the particle rapidity spectrum. 

The Fourier transform of the initial current can be easily obtained and they are related to j^{p^ ) of Eqs. (|78p and 
^ by 



i=l 

Similarly, the total Fourier component for Eq. is 

n 



(89) 



(90) 



where c{p^) is given by Eq. (|55p . As a consequence, Eq. (|64p for the momentum distribution of produced particles is 
modified to be 



^^tot ^ Ctot(p^)Ctot(p^) 

dpi 



2pO 

n 

E cos{p^(xj -X])} 



2pO 



(91) 
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The rapidity distribution for the fragmentation of n strings is 

dNtot 



dy 



^ 2CtCj cos{(m smhy){xt - Xj)} 



For the Fourier coefficient given by Eqs. ((501) . we obtain 



dy 



^ 2QCj cos{(m smhy ){x^ - Xj)} 



sinh^ ^ ' 



where ^ = i^^Tr^rnsinhy/Sv^ is given by Eq. (|83p . This equation can also be written as 



dNu 



dy 



^ 2CjCj cos{(m sinhy) {xi - X])} 



i>3 



dN^w 



dy 



where dNsingic/ dy = I'^^^/sinh^ f . The total energy of the n-string system is given by 



= J dym CiCj cosh y 



^ 2CjCj cos{(m sinhy) {xt - X])} 



sinh^ i ' 



which provides a relation between and cr. 

The Wigner function for the fragmentation of n strings is now modified to be 



/tot(a;, fc) 



47rfc0 



n 



iq 



2^ 



We can simplify the above expression by noting that 



giga;g-i(fc+(3/2)xi+i(fc-<?/2)xj _ giq{x-Xij)-iHXi-Xi} 

where Xij is the average of the initial longitudinal positions of strings i and j, 

X^i, = +Xj)/2- 

Therefore, we get 

n n n 

^ QQe'<l-e-'i'-+1mx.+^{k~<lmx, ^Y^e^'ii-x.) ^Y^QQ^^<l{--x.,)2QQcos{k{x^-X3)}■ 

i,j=l i i>j 



(92) 



(93) 



(94) 



(95) 



(96) 

(97) 
(98) 

(99) 



We can thus write the total Wigner function of the produced particles in the fragmentation of n strings as 

n n 

ftot{x, fc) = ^ /singic(a; - Xi,k) + ^ 2C,Q cos{fc(x, - Xj)} /single (2: - Xi],k). (100) 

i—l i>j 

The corresponding Wigner function ftot{x,y) can then be obtained as ftot{x,y) = ftot{x,k) evaluated at fc = 
m sinh(y). 

We note from these results that there are important interference effects between identical strings. The interference 
appears in the form of a cosine function C^iQ cos{fc(xi — Xj)} that depends on the momentum of the produced particle, 
the initial spatial separation between the strings, and the senses of the strings. The dependence on the string separation 
is similar to the interference in the emission of particles from different sources in intensity interferometry |87l |. 

In the case when the strings are far apart such that \xi — Xj\ is large, the cross terms in the Wigner function with 
i ^ j 'm Eq. (jlOOp oscillate rapidly about zero (for non-zero modes) . The cross terms will provide only small average 
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contributions and the direct terms dominate. In that case when strings are far apart, /tot(2;, k) w /singicl^; — Xii k) 
for |fc| 7^ 0, which is the sum of independent Wigner functions for separated strings. 

In the other extreme, all strings are located at the same spatial point. Among these n strings, there are n_|- strings 
with the C = 1 sense and n_ = n — n_|- strings with the opposite sense, C = ~1- Then, when all the strings occur at 
the same point, we have 

dNtot/dy - (n+ - 7i_fdN,,,,^iJdy, (101) 

and 

ftot{x, k OT y) = (n+ - ?i-)^/singio(a::, k or y). (102) 

(103) 

If all the strings are aligned in the same direction then dNtot/dy = v?' dN single I dy. On the other hand, if n_|_ = n_, 
then dNsingic/ dy = and ftot{x,y) = 0. In general, the greater the difference |n+ — the larger is dNtot/dy and 
ftot{x,y). The strength of the Wigner function of one string of one sense is canceled by the Wigner function of the 
string of the opposite sense. This cancellation of the string strength cannot be complete if the number of strings is 
odd. 



VIII. WIGNER FUNCTION IN THE FRAGMENTATION OF IDENTICAL STRINGS 




FIG. 3: (Color online) The Wigner function f{x,y) for two strings with the same sense, ^1^2 = 1- The strings are initially 
located at x — —0.1 and x = 0.1 fm indicated by the two arrows, with a separation of Ax = 0.2 fm. Fig. 3(a) is for t = and 
Fig. 3(b) is for t = 0.4 fm/c. 



As an example, it is instructive to evaluate the Wigner function for a few cases to see its dependence on the 
separation between the strings and their relative sense of directions. We show the time variation of the Wigner 
function for two strings each of which, if separated, has characteristics the same as those in Fig. 2. We examine first 
the case of two identical strings with the same sense so that C1C2 = 1- In Fig- 3, the origins of these two strings 
are initially located at xi ~ —0.1 and X2 ~ 0.1 fm with a separation Ax = 0.2 fm. In Fig. 4, they are located at 
Xi = —0.04 and X2 = 0.04 fm with a separation A.x = 0.08 fm. The locations of the origins of the strings are indicated 
by two thick arrows in Figs. 3 and 4. 

We note the following features from the dynamics of the Wigner function from Figs. 3 and 4. 
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FIG. 4: (Color online) Wigner function f{x, y) as a function of x for different values of y for two strings located sX x = —0.04 
and X = 0.04 fm, as indicated by the two arrows, with a separation of Ax = 0.08 fm. Fig. 4(a) is for t = G and Fig. 4(b) is for 
t = 0.4 fm/c. 



1. Similar to the case of a single string, all bosons of different momenta in different regions of the rapidity plateau 
are present at t = 0. This picture is in contrast to the momentum-space-time ordering of the classical picture 
of particle production. 

2. The magnitude and the width of Wigner function of the zero- mode with ?/ = are quite large when C,iC,2 = 1- 
The magnitude of the peak is about three times that of Fig. 2. For this case of two strings, there are three 
contributions, corresponding to two contributions centered at the origins of the strings, plus an additional 
contribution at a; = 0, as given by Eq. (jlOOp . For the y=0 zero mode, all these three pieces merge together to 
form a broad peak at a; = that is independent of time. 

3. For the low rapidity case with y = 1, the different contributions to the Wigner function merge into a single peak 
for both separations of Ax — 0.2 and 0.08 fm. 

4. For higher rapidities with y >2, the different contributions to the Wigner function merge into a single peak for 
the case of A = 0.08 fm, but splits into two peaks for Ax = 0.2 fm. This shows that the Wigner functions for 
many strings tend to merge together into a single structure, when the separation between strings is small. 

5. As a function of time, the peaks of the Wigner function move in the positive x direction for positive values of y. 
They move to the negative x direction for negative values of y. The low-momentum Wigner functions with y ~ 1 
remain to be merged as they propagate. For the case of y > 2 and Ax = 0.2 fm, there appear two propagating 
peaks corresponding to the propagation of the produced particles from the two different strings. For the case 
y > 2 and Ax ~ 0.08 fm, the Wigner functions from the two strings remain merged during their propagation. 

It is instructive to evaluate next the Wigner function for the case of two identical strings with the opposite sense 
so that ^1^2 = —1- We consider again two strings each of which, if separated, has characteristics same as those in 
Fig. 2. In Fig. 5, these two strings are initially located at xi = —0.1 and X2 = 0.1 fm with a separation Ace = 0.2 fm. 
In Fig. 6, they are located at xi = —0.04 and X2 = 0.04 fm, with a separation Ax = 0.08 fm. The positions of the 
origins of the strings are indicated by the two thick arrows in Figs. 5 and 6. 

The results in Figs. 5 and 6 show that in the fragmentation of two strings with opposite senses, the magnitude 
of the Wigner function is much reduced from the fragmentation of two strings with the same sense. The zero-mode 
Wigner function is negative at x = 0, in contrast to the positive definite property of the Wigner function for a single 
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FIG. 5: (Color online) The Wigner function f{x,y) as a function of x for different values of y for two strings with opposite 
senses, 1^11^2 = —1- The strings are initially located at x = —0.1 and x — 0.1 fm indicated by the two arrows, with a separation 
of Ax = 0.2 fm. Fig. 5(a) is for t = and Fig. 5(b) is for t = 0.4 fm/c. 



string or for two strings with the same sense. There are oscillations of the Wigner function for different momenta at 
t = 0. For i > 0, the dominant peaks of the Wigner function propagate according to their momenta. 

As the separation between the strings becomes smaller, the Wigner function of opposite senses tend to cancel and 
the magnitude of the Wigner function decreases, as shown in Fig. 6. As we remarked earlier, the Wigner function is 
reduced to zero in the limit when the two identical strings with opposite senses coincide. 



IX. DYNAMICS OF THE WIGNER FUNCTION IN THE FRAGMENTATION OF INDEPENDENT 

STRINGS 



The last two sections deal with the fragmentation of identical strings for which interference effects are present. In 
a nucleus-nucleus collision, the number of nucleon-nucleon collisions is large in a local neighborhood. We expect that 
each nucleon-nucleon collision produces at least two strings spanned between the quark of one nucleon and the diquark 
of the other nucleon. The types of strings that are formed are however very large in numbers as each string will be 
characterized by the color and the flavor of the quark and the diquark at its two ends, and by the color gauge field 
component a of A'^'^ the string can excite. As a consequence, the probability of a neighboring string to be identical 
is of order \/[N'^N f{q)N f{qc[){N'l — 1)] which is quite small. The production process in string fragmentation in a 
nucleus- nucleus collision is likely to involve the fragmentation of independent non- identical strings. It is appropriate 
to examine the space-time dynamics of produced particles in the fragmentation of many independent strings. 

In a nucleus- nucleus collision in the ccntcr-of-mass system, the average longitudinal separation between neighboring 
strings is (i/7 where d ~ 1.9 fm is the average nucleon-nucleon separation for a nucleus at rest and 7 = ^s„„/mjv is 
the Lorentz contraction factor. On the other hand, the width of the Wigner function for the zero mode with y^O is 
of order h/rriT^T, as indicated by Eq. (|86p . Whether the Wigner function will be separated or will overlap depend on 
whether d/j is greater or less than h/rriT^T- 

a d/j > h/niT^T, strings are separated. (104) 
If '^/t f^/fn-rrT, strings overlap. 

As the pion transverse mass mTrT is of order 0.5 GeV, the Wigner functions will be separated from each other if the 
collision energy per nucleon in the ccntcr-of-mass system is y/s^^ ^ 5 GeV and will overlap each other if ^s„„ ^ 5 
GeV. 
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FIG. 6: (Color online) The Wigner function f{x,y) as a function of x for different values of y for two strings with opposite 
senses, (1(2 = —1- The strings are initially located at a; = —0.04 and x = 0.04 fm indicated by the two arrows, with a separation 
of Aa: = 0.08 fm. Fig. 8(a) is for t = and Fig. 8(b) is for t = 0.4 fm/c. 




FIG. 7: (Color online) Schematic space-time evolution of the Wigner function f{x,y) of produced particles from many 
independent strings. Fig. 7(a) is for t = Q, 7(b) and 7(c) are the schematic configurations at subsequent times. 
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For each string, particles of different rapidities in different regions of the rapidity plateau are produced at the moment 
of string fragmentation. They propagate toward opposite longitudinal directions according to their rapidities. We 
can depict the dynamics of the density of Wigner functions as those of one of the rectangular regions in Fig. 7. 

How the partons evolve after production in a heavy-ion collision at high energies is a complex problem which will 
require extensive future investigations. Our experience with the phase space dynamics of nucleons in nucleus- nucleus 
collisions as investigated earlier in (ssl . Isoj may furnish an approximate guide to speculate on the behavior of partons 
in high energy nuclear collisions. 

We can discuss first the case of y/s^ ^ 5 GeV in the fragmentation of independent and separated strings as 
shown in Fig. 7. Produced particles in these separated strings will stream from one point to another point as time 
proceeds. Their streaming will lead to three important effects. First, they will replenish those particles that have 
just streamed away from the point of production, as is depicted in Fig. 7(b) and 7(c). Secondly, as a result of 
their streaming, they will encounter at the same longitudinal point produced particles streaming from neighboring 
points in the opposite direction. Such encountering will lead to collisions. The discussion of the collisions of the 
produced particles necessitates the inclusion of the transverse degrees of freedom. It is necessary to restore the 
transverse momentum distribution of the produced particles by using information on the transverse eigenfunctions 
of the quarks and antiquarks. These collision of produced particle will convert longitudinal momenta to transverse 
momenta, leading to medium particles with smaller longitudinal momenta but greater transverse momenta in the 
approach to equilibrium. Thirdly, the produced particles will likely be subject to strong residual interactions after 
production. The interaction is stronger, the more the strings are closer packed together at higher collision energies. 
Fourthly, because of their mutual interaction and the nature of color confinement, those produced particles reaching 
the boundary region will form a surface in the boundary region, against which produced particle arriving later will 
collide. Part of these particles will continue to propagate forward as a result of the collision with the longitudinal 
boundary and will stretch the surface to a greater longitudinal extension. However, part of the colliding produced 
particles will be reflected from the surface and will stream in the backward direction, replenishing some of the partons 
that have streamed backward, similar to the description in Fig. 1 of (soj . in the approach to equilibrium. 

We consider next the heavy-ion collisions at very high energies. For a collision energy of 100 GeV per nucleon 
at RHIC, the average longitudinal separation between strings is of order 0.019 fm which is much smaller than the 
width of the Wigner function, h/niT^T- The nucleus-nucleus collision leads to a large overlap of the fragmenting strings. 
Because of the strong overlap of the strings, the produced quanta will propagate to the longitudinal boundary together 
in nearly the same speed and nearly in unison. Because of the strong overlap and the Lorentz contraction, there is 
a large density and large interaction among produced particles after they are produced. The produced particles 
will propagate according to their longitudinal rapidities. Partons with large rapidities will reach the boundary first, 
forming a surface region. The interacting partons near the surface will collide with the confining and moving surface. 
As a result of the collision with the boundary, part of these partons will continue forward, and they will stretch the 
surface boundary to a greater longitudinal extension to keep it a continuously moving boundary. Because of the 
confining interactions at the moving boundary, some of the other partons will come to the classical turning point at 
the moving boundary. They will then be refiected from the surface and will flow in the backward direction. The 
reflection of confining partons leads to two important effects. First, the reflected partons will replenish some of 
the partons that have streamed backward, as depicted in Fig. 1 of Ref. [s^. Secondly, the reflected partons will 
encounter partons streaming in the opposite direction. Such encountering will lead to collisions which will convert 
longitudinal momentum to transverse momentum, leading to medium particles with a smaller average longitudinal 
momentum. The collisions may relax the rapidity distribution from an initial rapidity plateau [3l] | to a Gaussian 
rapidity distribution at the end of nucleus-nucleus collision [o^, [Ml, [111 . Whether or not such a scenario indeed occurs 
will require more future analyses. 

X. DISCUSSION AND CONCLUSIONS 

We investigate the space-time behavior of produced particles in the fragmentation of a string in QED2. We obtain 
a relation between the rapidity distribution and the Fourier transforms of the initial fermion current. We find that the 
rapidity plateau at high energies arises as a consequence of a localized initial fermion current. The Wigner function 
of the produced particles show that particles with momenta in different regions of the rapidity plateau are present 
at the moment of string fragmentation, with the width of their spatial distributions decreasing with an increasing 
center-of-mass energy. The Wigner function exhibits the effects of interference in the fragmentation of many identical 
strings. 

Because QED2 mimics many features of QCD4, it is useful to examine the circumstances in which a QCD4 sys- 
tems can be approximated by QED2. We show first how QCD4 with transverse confinement can be approximately 
compactified as QCD2 with a transverse quark mass obtained by solving a set of coupled transverse eigenvalue 
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equations. Furthermore, in the hmit of the strong couphng and the large number of flavors Nf, QCD2 admits 
Schwinger QED2-type solutions. It is for these reasons that QED2 can mimic many important features of QCD4, 
including the properties of the proper high-energy rapidity plateau behavior, quark confinement, charge screening, 
and chiral symmetry breaking. In the absence of rigorous non-perturbative QCD4 solutions, it is therefore reasonable 
to study the particle production process phenomenologically using Schwinger's QED2 model, with the boson quanta 
of QED2 considered as analogous to the boson quanta in QCD. 

There are however important differences which must be kept in mind in the discussion of the dynamics subsequent 
to the fragmentation of strings. Free bosons are the quanta of QED2 while strongly interacting gluons and quarks 
and relatively weakly interacting hadrons are quanta of QCD4 at different temperatures of the QCD system. In a 
nucleus-nucleus collision at high energies, the particles produced after the fragmentation the strings will be subject to 
strong interactions with the dense medium of produced particles in their vicinity after production. These interactions 
are important as the strings overlap strongly owing to the Lorentz contraction. The large number of colors and flavors 
make it likely that the overlapping strings are of different types and do not interfere. As a consequence, the density 
of the produced particles accumulatively increases as the density of strings increases. It will be necessary to include 
these strong interactions between the produced particles in order to describe better the subsequent dynamics. It will 
also be necessary to restore the transverse degrees of freedom for a better description by using information on the 
transverse eigcnfunctions of the quarks and antiquarks and their transverse momentum distributions. Our formulation 
of the approximate compactification between transversely-confined QCD4 into QCD2 facilitates such a restoration. 

The main feature inferred from our present analysis is that in a string fragmentation, boson field quanta with 
momenta in different regions of the rapidity plateau are present in the initial Wigner function at the moment of string 
fragmentation. This peculiar behavior arises from the quantum effects of the vacuum structure. Particles in the 
vacuum of the Dirac sea participate in the interaction as the valence fermions, antifermions, and gauge fields change 
with space and time. The non-perturbative self-consistent response of particles in the Dirac sea leads to the production 
of boson quanta with momenta in different regions of the rapidity plateau at the moment of string fragmentation, in 
contrast to the classical picture of string fragmentation where there is a momentum-space-time ordering of produced 
particles. 

Our analysis of the dynamics of the Wigner function places it in the class of initial-state-interaction description 
of quantum chromodynamics. Our result that partons in different regions of rapidities over the rapidity plateau are 
produced at the moment of collision is in line with many other well-justified initial-state-interaction models, such 
as the parton model [6l|, the Drell-Yan process [g^, the multi-peripheral model [6^, and the dipole approach of 
photon-induced strong interactions [6^ . These are models in which the constituent particles or produced particles 
are present at or before the collision process at t = 0. 

With regard to the momentum kick model which motivated the present analysis, the above results may resolve one 
of the puzzles concerning the possible occurrence of the particles with large longitudinal momenta in the early stage of 
the nucleus-nucleus collisions. Partons of different rapidities are present in the initial parton momentum distribution 
at fragmentation. The jet produced in one of the nucleon-nucleon collisions can find partons of large rapidities in 
the early environment. There can be collisions between the jet and these high-rapidity partons, which show up as 
ridge particles in coincidence with the jet. Just as the presence of antiquark partons from the quark-antiquark sea 
is revealed by the occurrence of the Drell-Yan process in the initial-state-interaction picture, so it is here that the 
presence of the rapidity plateau in the early parton momentum distribution is revealed by the occurrence of large 
rapidity associated particles in coincidence with the near-side jet in the PHOBOS experiments. 
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